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3 The Moon's changeable aspect during a lunar eclipse is largely attributable 

4 to variations in the refracted unscattered sunlight absorbed by the terres- 

5 trial atmosphere that occur as the satellite crosses the Earth's shadow. The 

6 contribution to the Moon's aspect from sunhght scattered at the Earth's ter- 

7 minator is generally deemed minor. However, our analysis of a published spec- 

8 trum of the 16 August 2008 lunar eclipse shows that diffuse sunhght is a ma- 

9 jor component of the measured spectrum at wavelengths shorter than 600 

10 nm. The conclusion is supported by two distinct features, namely the spec- 

11 trum's tail at short wavelengths and the unequal absorption by an oxygen 

12 coUisional complex at two nearby bands. Our findings are consistent with 

13 the presence of the volcanic cloud reported at high northern latitudes fol- 

14 lowing the 7-8 August 2008 eruption in Alaska of the Kasatochi volcano. The 

15 cloud both attenuates the unscattered sunlight and enhances moderately the 

16 scattered component, thus modifying the contrast between the two contri- 

17 butions. 
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1. Introduction 

18 The classical theory of lunar eclipses accounts for refraction, differential absorption and 

19 focusing to explain the Moon's aspect during an eclipse [Link, 1962]. Link's classical the- 

20 ory has been subsequently perfected and used to investigate the composition and aerosol 

21 loading of the Earth's atmosphere [e.g., Ugolnikov and Maslov, 2008]. Aerosols play a 

22 critical role in the interpretation of lunar eclipses as their content, distribution and opti- 

23 cal properties are largely unpredictable. Volcanic eruptions and meteor showers are two 

24 natural sources of aerosols with the potential for perturbing the atmosphere and, in turn, 

25 the aspect of the eclipsed Moon [Keen, 1983; Vollmer and Gedzelman, 2008]. Occasionally, 

26 large wildfires may also perturb the atmosphere [Fromm et ai, 2010]. 

27 Garcia Munoz and Palle [2011] have revisited the lunar eclipse theory and estimated the 

28 impact of volcanic aerosols on the spectrum of sunlight at the eclipsed Moon. Aerosols may 

29 substantially attenuate the direct sunlight while simultaneously enhancing somewhat the 

30 scattered contribution. The latter depends strongly on the capacity of aerosols for forward- 

31 scattering the incident light and, consequently, on the aerosols' size. The spectroscopic 

32 characterization of the sunlight reflected from the eclipsed Moon takes the investigation 

33 of lunar eclipses farther than allowed for by photometry and the traditional color indices. 

The 7-8 August 2008 eruption of the Kasatochi volcano (52.17°N, 175.51°W, Aleutian 

35 Islands, Alaska) ended a period of global low stratospheric aerosol amounts. The three 

36 main explosions recorded over two days plus the release of gas that followed for hours 

37 delivered into the atmosphere ~1.5 Tg of SO2 [Waythomas et al, 2010], which is ~30 

38 times less than the SO2 injected by Pinatubo in 1991. The plume of gas and ash rose up 
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39 to ~14-18 km and drifted eastward carried by jet winds, spreading rapidly over North 

40 America, Greenland, and the North Atlantic Ocean. The cloud was spotted above Europe 

41 on 15 August, one week after the eruption [Martinsson et ai, 2009]. 

42 Our paper shows that a published spectrum of the Moon in umbra during the August 

43 2008 lunar eclipse contains sunlight scattered at the Earth's terminator. We argue that the 

44 Kasatochi eruption is the most plausible origin for the abnormally elevated atmospheric 

45 opacity needed to explain the observation. Vidal-Madjar et al. [2010] have published a 

46 spectrum of the August 2008 lunar eclipse, but covering only the penumbra. 

2. Data 

47 We use data of the 16 August 2008 lunar eclipse obtained with the ALFOSC instrument 

48 mounted on the Nordic Optical Telescope at the Observatorio del Roque de los Muchachos 

49 (La Palma, Spain) and presented by Palle et al. [2009]. The dataset comprises spectra 

50 of the Moon in umbra (21:36UT), penumbra (22:11UT), and out of echpse (23:09UT). 

51 Palle et al. [2009] derived a lunar eclipse spectrum from the ratio of umbra and penumbra 

52 spectra. The ratio cancels out the solar spectrum and the telluric signature of the Moon- 

53 to-telescope optical path. What remains is the imprint of the limb-viewed terrestrial 

54 atmosphere (averaged in a particular way over the terminator) on the sunlight that reaches 

55 the Moon in umbra. Our analysis sets out from the 400-900 nm published spectrum. 

56 The solar elevation angle, e, is the geocentric angle between the incident sunbeam 

57 direction and the direction from the Earth's centre to the lunar disk parcel targeted by 

58 the telescope. We have that e~0.34° for the slit projected on the Moon. The structure of 

59 an umbra spectrum is very sensitive to e [Garcia Mufioz and Palle, 2011]. 
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3. Evidence of scattered sunlight 

3.1. The short- wavelength tail of the spectrum 

60 Model predictions for e~0.34° and a broad range of aerosol loadings show that the 

61 spectrum of sunlight directly transmitted through the atmosphere is typically 2-3 orders 

62 of magnitude fainter at 400 nm than at 600 nm [Garcia Munoz and Palle, 2011]. This 

63 is at odds with the eclipse data, which show that the measured spectrum is roughly flat 

64 to within a factor of 2 shortwards of 600 nm and that the fluxes at 400 and 880 nm are 

65 in a ratio of ~1:20. It thus means that direct sunlight is not the only contributor to the 

66 measured spectrum. Garcia Munoz and Palle [2011] note that a flat spectrum at short 

67 wavelengths indicates that diffuse sunlight dominates locally over direct sunlight. 

3.2. The (02)2 bands at 577 and 630 nm 

68 We fltted synthetic curves of the form niexp(— r^) to the measured spectrum from 

69 550 to 660 nm. The curves include absorption by H2O, O3, O2 and the (02)2 coUisional 

70 complex. One term, Tcont=Z]fc=o Cfe(A*/A)'^, with 600 nm, accounts for a continuum 

71 baseline. Thus, each curve contains up to ten degrees of freedom, namely, flve c^'s, 

72 integrated columns for H2O, O3 and O2, and, optionally, one integrated column for each 

73 of the [X^E-{0)]2^a^Ag{0)WAg{l) and [X^Eg {0)]2^[a^ Ag{0)]2 bands of (02)2 that 

74 occur at 577 and 630 nm, respectively. For the temperature-dependent gas properties, 

75 the temperature was flxed at 225 K. The synthetic curves were properly degraded and 

76 resampled. The minimization of X^=J2j{^ ~ -^fit(Aj)//exp(Aj))^, where /fit and /gxp are the 

77 synthetic and observed data, outputs the best flt parameters. 
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78 Figure (1) summarizes the best fits obtained from tliree separate strategies, eacli of tliem 

79 treating the 577 and 630 nm bands of (02)2 in a different manner. Strategy A fits the 

80 spectrum with null amounts of (02)2; B includes (02)2 and assumes the same integrated 

81 column for the two bands; and C allows for separate integrated columns for each of the 

82 577 and 630 nm bands. In the top panel, the solid black curves represent the measured 

83 spectrum, shifted in the vertical for comparison with the synthetic curves. The red solid 

84 curves are the respective A, B and C best fits. The bottom panel displays the residuals. 

85 Including (02)2 reduces notably the fit residuals. The fit improves further if the integrated 

86 column at 630 nm is about twice the column at 577 nm. The latter conclusion is the core 

87 of the second argument that proves the significance of diffuse sunlight in the measured 

88 lunar eclipse spectrum. Some comments on the robustness of the fitting procedure can be 

89 found in the Supplementary Material. 

90 Taking C as the optimal strategy, the conclusion is that average sunlight photons at 577 

91 and 630 nm follow different paths in the atmosphere. The direct trajectories of sunlight 

92 rays are dictated by the atmospheric refractive index, which does not change appreciably 

93 within such a narrow spectral interval. The amount of sunlight directly transmitted does 

94 however vary sharply with wavelength. Direct sunlight is more attenuated at 577 nm 

95 than at 630 nm due to the ~A~^ behaviour of the Rayleigh cross section and the closer 

96 proximity of the 577 nm band to the absorption peak of the O3 Chappuis band. 

97 We thus have to invoke sunlight scattered at the Earth's terminator to explain the 

98 measured spectrum. Garcia Munoz and Palle [2011] show that in a lunar eclipse the 

99 bulk of diffuse sunlight near 600 nm originates from above 15 km. In the stratosphere. 
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100 the (02)2 density, which drops with a scale height half that of the background density, 

101 is negligibly small. Foreseeably, the signature of the (02)2 bands in the diffuse sunhght 

102 spectrum is weak. 

4. Analysis and discussion 

103 Next, we generate model lunar eclipse spectra that reproduce the measured spectrum 

104 if a few reasonable assumptions on the loading and properties of airborne aerosols are 

105 introduced. The spectra, generated with the model described by Garcia Munoz and Palle 

106 [2011], contain both direct and diffuse components. Further details on the underlying 

107 model assumptions can be found in the Supplementary Material. 

108 The tracing of the direct sunbeam that reached the parcel of the Moon's disk tracked 

109 by the telescope reveals that the sunbeam intercepted the volcanic cloud formed in the 
no Kasatochi eruption, as seen in Fig. (2). It is expected that the direct sunlight component 

111 is more strongly affected by the volcanic cloud than the diffuse one, which originates from 

112 all terminator locations. This distinction is accounted for in the generation of the model 

113 spectra by assuming separate atmospheres for the calculation of each component. 

114 For simplicity, the model spectra are allowed only four adjustable parameters. These 

115 are /yg, a' and /03 for the direct sunlight component, and Tefr for the diffuse one. In the 

116 former, /^^ scales the reference aerosol extinction profile at 1.02 //m, a' is the Angstrom 

117 exponent to extrapolate the 1.02-/xm extinction profile to shorter wavelenghts, and /03 

118 scales the reference ozone profile. In the calculation of the diffuse sunlight component, r^G 

119 stands for a mean effective radius for aerosols at the terminator. The sulfate droplets of 

120 background aerosols in the quiescent atmosphere have reflf~0. 1-0.2 //m, whereas volcanic 
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121 ash particles with residence times longer than a few days may have refr's of a few microns 

122 [Bauman et ai, 2003; Munoz et al, 2004]. Sioris et al. [2010] report Tes's of ~0.6 /xm 

123 for September 2008, which are indicative of a perturbed atmosphere. It is unclear what 

124 the mean at the terminator was one week after the eruption. Thus, we explored a set 

125 of Teff from 0.1 to 2 /xm to bracket possible sizes. A large Tefr results in phase functions 

126 strongly peaked in the forward direction. 

127 For each res, one diffuse sunlight spectrum was produced. For each diffuse spectrum an 

128 algorithm seeks the /^q, a' and /03 values producing the best fit of the direct + diffuse 

129 model spectra to the continuum of the measured spectrum. The algorithm forces the 

130 (flux-uncalibrated) measured spectrum to match the model spectra at 875 nm. Figure (3) 

131 shows the best fit for refr=0.5 /xm and the values inferred for the other three adjustable 

132 parameters. The a parameter is the multiplicative factor to pass from the normalization in 

133 the graph to the Earth-to-Sun ratio as discussed by Garcia Munoz and Palle [2011]. Figure 

134 (I) in the Supplementary Material shows the best fits for the full r^s set. It is apparent 

135 the good a posteriori match of the O2 bands in all cases. The differences between the 

136 measured spectrum and the best fits are of a few percent longwards of 600 nm, but of 

137 ~50% near 500 nm. This is a consequence of fitting the measured spectrum with models 

138 that contain a reduced number of adjustable parameters. The residuals longwards of 700 

139 nm are mainly due to a known instrumental issue of uncorrected fringing. 

140 The /-yo values inferred point to heavy aerosol loadings with peak extinctions of ~10~^ 

141 km~^ in the atmosphere intercepted by the direct sunbeam. Comparable extinctions 

142 were reported on global scales for a few months after the Pinatubo eruption [Bauman 
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143 et ai, 2003]. Exponents a'~0-0.15 are indicative of large-size particles being carried in 

144 the volcanic cloud. The conversion from SO2 to sulfate droplets has an e-folding time 

145 of 20-50 days [Kristiansen et ai, 2010]. It is unlikely that one week after the eruption 

146 is enough time for large sulfate droplets to form. Thus, the a' values inferred suggest 

147 that the volcanic cloud carried sizeable amounts of unsedimented ash. Sioris et al. [2010] 

148 report small Angstrom exponents of ~0.5 in early September 2008. 

149 The inset of Fig. (3) shows the two component spectra near 600 nm. The diffuse 

150 spectrum is roughly flat and shows no evidence of (02)2 absorption. When the direct and 

151 diffuse model spectra are added, the 577 nm band becomes more diluted than the 630 

152 nm band, which translates into an effective integrated column at 630 nm larger than at 

153 577 nm. For the case in Fig. (3) the ratio is ~1:1.4, somewhat smaller than the ~1:2 

154 ratio inferred from the measured spectrum. One may generally state that comparable 

155 amounts of direct and scattered sunlight near 600 nm lead to larger (02)2 columns at the 

156 longer-wavelength band. 

157 Figure (I) in the Supplementary Material proves that good fits to the measured spectrum 

158 are possible for the full Ves set. This means that the measured spectrum accepts one 

159 quantitative interpretation for each rcff. In qualitative terms, though, the picture that 

160 we obtain is fairly consistent and indicates that the direct sunbeam was substantially 

161 attenuated by the volcanic cloud, which leads to an enhanced contrast of the diffuse 

162 component. For future efforts, we suggest that the flux calibration of the undivided 

163 spectra might help break the degeneracy. 
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164 A comment is to be made regarding the Ring effect and the structure seen in the mea- 

165 sured spectrum shortwards of 540 nm. The Ring effect refers to the smearing of solar 

166 Fraunhofer hues that occurs in the spectrum of sunhght scattered in the atmosphere 

167 [Grainger and Ring, 1962]. In the Earth's atmosphere, the Ring effect is due to rotational 

168 Raman scattering by N2 and O2 [Kattawar et ai, 1981]. Raman scattering redistributes in 

169 wavelength part of the incident photons, the redistribution being more evident where the 

170 incident solar spectrum shows the sharpest lines. The ratio of scattered to unscattered 

171 sunhght spectra reveals the Ring effect as a fiUing-in of the solar hne cores. The detection 

172 of the Ring effect in the eclipse data would mean a further confirmation of scattered sun- 

173 light. The measured lunar eclipse spectrum shows that ripples do occur in the Fraunhofer 

174 region. The sign of the structures is however inverted with respect to what the Ring effect 

175 would produce. The inspection of the undivided umbra spectrum shows that the solar 

176 lines are unexpectedly deep, probably due to the limitation in the subtraction of the sky 

177 spectrum at these wavelengths, where the signal-to-noise ratio is the lowest. Thus, the 

178 structure seen in the measured spectrum cannot be attributed to the Ring effect. Further, 

179 a few quantitative arguments allow us to deem as minor the impact of the Ring effect on 

180 the measured lunar eclipse spectrum. Following Kattawar et ai [1981], the filling- in for 

181 forward-scattered sunlight is ^^2.5% of the continuum Rayleigh-scattered by the gas. In 

182 the conditions explored here the filling-in would be undetectably small because the sun- 

183 hght scattered by the gas contributes less than a few percent to the net sunhght scattered 

184 by gas and aerosols together. 



DRAFT 



June 16, 2011, 12:29ain 



DRAFT 



GARCIA MUNOZ ET AL.: VOLCANIC AEROSOLS IN A LUNAR ECLIPSE X - 11 

185 Pyro cumulonimbus (pyroCbs) is a recently-coined term to designate convective activity 

186 triggered or sustained by wildfires [Fromm et ai, 2010]. In extreme events, pyroCbs inject 

187 smoke and biomass-burning particles into the troposphere and lower stratosphere and alter 

188 the global aerosol loading. PyroCbs may result in aerosol extinctions ~10~^-10~^ km~^ 

189 well above the tropopause, opacities that are often associated with volcanic clouds. It 

190 would be difficult to differentiate the impact on the eclipsed Moon of one such event from 

191 that of a volcanic eruption. To our knowledge, no extreme pyroCbs were reported in 

192 the weeks preceding the eclipse, a period that was monitored with unprecedented detail. 

193 Thus, if any, the contribution in the eclipse of pyroCbs blended with that of the Kasatochi 

194 cloud. 

195 The Perseids is one of the most copious meteor showers, running yearly from late July 

196 to late August. In 2008, its peak of activity occurred near 13 August. Despite recent 

197 work [Mateshvili et ai, 1999; Renard et ai, 2010], there are significant uncertainties on 

198 the optical properties of the atmosphere perturbed by meteor showers. Mateshvili et al. 

199 [1999] report two-fold enhancements with respect to pre-shower values in the twilight 

200 brightness above 20 km during the Leonids in 1998. Assuming that both meteor showers 

201 are comparable and that the brightness enhancement translates into a similar increase in 

202 stratospheric opacity, the effect of extraterrestrial dust would be more than one order of 

203 magnitude smaller than that by the volcanic perturbation. Thus, the effect of meteoroid 

204 dust in the measured spectrum is likely masked by the volcanic perturbation. 

205 We have shown that the lunar eclipse spectrum published by Palle et al. [2009] was 

206 affected by sunhght scattered at the Earth's terminator. We offered theoretical arguments 
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207 that hint at the Kasatochi eruption as the most plausible origin for the atmospheric 

208 perturbation needed to explain the observations. Future observations will allow us to 

209 compare lunar eclipse spectra obtained in different atmospheric conditions. In a broader 

210 context, it is worth mentioning that the retrieval of globally- averaged atmospheric optical 

211 properties is a relevant exercise towards the future characterization of transiting Earth- 

212 hke extrasolar planets. As a corollary, we may state that the color of the lunar disk in 

213 umbra during the 16 August 2008 lunar eclipse was partly caused by diffuse sunlight. 
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216 for constructive comments. 
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Figure 1. Top: Best fits (red) to the measured spectrum (black) from 550 to 660 nm. Tlie 
dotted green line are the best fits divided by exp (— T(02)2)- The comparison of the dotted and 
solid curves makes explicit the contributions from the coincidental in position, albeit distinct 
in nature, O2 X(0)^6(2) and (62)2 X(0)+X(0)-)-a(0)+a(0) absorption bands near 630 nm. 
Bottom: Fit residuals. For B and C, the dotted red curves are the (©2)2 contributions. In C, we 
infer an optimal ratio for the 577:630 nm integrated columns of ~1:2. 
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Figure 2. Solid: Projected mid-section trajectory of the sunbeam that reaches the lunar disk 
targeted by the telescope at 21:36UT on 16 August 2008. Overplotted, the SO2 cloud (a usual 
volcanic cloud tracer) on 15, 16 and 17 August (red, green and blue, respectively) according to 
AURA/OMI data (downloaded from the Giovanni online data system, developed and maintained 
by the NASA GES DISG). The sunbeam's closest approach to the Earth's surface occurs in the 
North Atlantic region. The local tropopause is at ~10 km. Dashed: Line of the terminator. 
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Figure 3. Best model fit for refr=0.5 fim (red) to the measured spectrum (black). The model 
spectrum contains contributions from direct sunlight (green) and diffuse sunlight (blue). The 
inset is a zoom of the region near 600 nm. The algorithm aims the fit of the continuum away 
from O2 and H2O bands. The H2O bands were fitted separately after the fit to the continuum. 
Figure (I) in the Supplementary Material shows the best fits for the full set of Tcs values. 
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Supplementary material 
AdditiioiLai cmnmenls aa Out robisstmfss eitbs; StihsgpFoceduie. af Seetkm^ 

Our fitting algcriflmmsK tie laboratoii' nseaajiEiDent^ for (0;)j: cros sections ly Gieieiiblaft et al. 
(1990). Later espeiiiujeuta [K'ewiLhaiu and Ballard. Kms and Ubachs. 1999] and transmitted 
Quooliglit uQeasaremeuts {Wkgoer et al., 2002] haw confimied that tbe Oieeiiblatt al. (1990^ 
cross -sectioiis are couect to within a few percent Ufe verified that using the Nemahana and BalSard 
{1998) cross sections does nai naodi^'' our conclusion that Etiateg\- C peiforais better than stnate^s. 
AandB. The condusiau lemain^ also umcIuiiDged if auaEteanative teuopeiatiiiemtiieiange 200-2^ 
K i£ used to ei'aluaile the teoqieiatuie-dependait optical poperties of tbe iat^oitiing gases. 

Additjonal mforioatiion on the spectra simiiladDiis of Sfctiion A 

A lijoiDU^ desciiption of fiie model nsed in Section 4 can be found in Garcia Midjoz and Palle 
{2011). \%iticallY-iesoived optical piopaties of the atmosptsffe are mput to the model The model 
al&o requiies the scattering phase fiinctions of all aiibome particles . Then, for a selected solai 
eEevadon angle, it produce both, direct and diffuse spectia of the sirnhght am^'ing at the lAmn 
noinaaliaedto the net solai inadianceL 

We follow Garcia ^Amoz and PaUe (2011) for the taicperafMre anddiensities of the spectRJscopicalty 
acti\'e naolecules in the spectral region of innestiEHtion (maiiilv HiO, O: and O.h). The composite 
gjcfcal map of cloud top heidts retrieved from GO\5E-2 data (ICSU Worid Data Centa^ for Remote 
Sensing of the Atmosphere. http:j"wdc..dlr.deO shows for 16 Augn^ 2003 clouds with top heidits 
between 4 and 10 km at 50-70"N latitudes. In the direct sunlight calculation, we assume that clouds 
block the Fmnliriit below 6 koi A precise choice of the cloud top height is not critical because for 
elevated aerowl loadings the atmosphae at that altitude is opticalty thick in limb ^ihiing. Most of 
the diSuse sunlight is scattered fi'om hiche~ altitudes, and the choice of the cloud top hed^t in the 
calculation of the dijESuse coDuponsit is less iiqjcatanl 

Several woriis haw imJiKtigated ihe optical pnopeities of the post-enqrtion atmospheie [e.g., 
Bourassa et al.. 2010; HofiBuann et al., 2010]. Sioris et al. C2010) ha^-e paroduced 1.02-}mi extinction 
profiles fi^om so'lar oceultation data obtained with the Atmospheiic Chemis^' ExpeiuQent (ACE) 
Imager. The September 200S ACE Imager estinctioii profile peaks at altitudes of ^-9 and 15 tm, and 
has an integrated optical thickness above 6 km of -X} .02. notably laigs than the optical thitiness for 
liK same month firan 2004 to 2007. It must be noted that tte Septenira: 2O0® ACE Imaga" 
estinction profile Chereafter refeired to as y,, Qj^ z)) wasi inferred firom data gathered 2-7 is'edcs after 
the eclqjse and includes measurements taken o^'Cr i»e entire northem hemisphere. ThiuE, -ft, z) is 
not dir«£y applicable to flie inteipFetation of the mea&iued kmar eclipse spectnnu. 

Ear BUT direct sunb^ calculalion5, we take as ortinction profile by asosols y (>, 2}=^ ^, 70 (Jo, 2) 
(X/V)^, wheie f ^1 is a mere scaling &ctor and q' is a so-called An^Hom exponent. Typically, a' 
ranges from 0, fcs: parficles much larger than the inddenl phottsi wawlensfth, to 4, in the limit of 
molecule-sized particles. IJfe also allow for ozone extinction in the 400-^50 nm leson of the 
Cliappiiis band. For this, we take the referaice ozone profile of Garcia Munoz and Palle {2011) and 
scale it by an additional parmieter f,:^. A battery of 17x53i>^25 direct simhght spectra was created 
to saK5ile fiie f - a' - tni parameter -space in the [0, 8I-[0. 2]-{0. 2] interi-als. Tbese spectra are 
mainly intoided to reproduce the continuum of the meastued spectnnu awi>' fiom mdecular bands. 

As a disck, it is worth noting that the \'alue3 lElriei^ for foj in the fits of the Section 4 range fl"0(m 
0.E5 to 0.95, which mean {for a reference piofile of -360 Dobson UnitE of ^nerticaDy-inle^ted 
OEone) ozone columiE bel?,™! 306 and 342 DU. The total ozone nsap for IS August 200S {World 
OziHie and Ultraviolet RadiatioQ Data Centre, WOLTK, hlIp:'>'ft'w\v.woiidc.Qrg) shows that in the 
North Atlantic region near Greeailand, where the simbeam intscepts the terminator, the ozone 
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column is 300-350 DU. Thus, the ozcmje colmtiiis lepoiled at WOUDC and ofur own letrievied 
cohmms ire in reasonable agrmuent 

Our diffuse simli^ calculations use Hie September 2007 ACE Imager exdnctian pdofile as 
rqreesailative of background aerosols. The paofile is sonoeigdiat arbitrarily muftiplied hy 1 .5 to 
accoimt for some post-eruption siobfil aerosol enhancement. Shoitvrards of 1.02 fim. we infer the 
estinctioiiE by means of an Angstrdni-t}'pe law with an exponent of 1 ^jpropiiate to a weakly-to- 
moderately perturbed atmospha'e [Haji'ashida and Hbiikawa. 2001] . The total ozone map for 16 
August 200S shows that at ■flie Earlii's taminator the awiage ozsne column vi-as less than in oui 
refcienoe profile. Thus, for the diflflise sunh^ calculations we ha^-e le-scaled the referenoe profile 
to a total ozone cohmm of -250 Dobaon LMts. 

The phase function of autrome particles at the small scattsing angles that occur for tbrwirdr-, 
single-scittered sunlight photons is critical in ibs diSiise probleuL MiUipiLe processes, which 
include sedimentatiorL ack'eciion, condaisation and coagulation, go\™i the aerosols' ewlving size 
distriibudon in the weeks and montiLS aAer m eruption. It is not £iowii a prtoh the average global 
value of r^F on the date of the edipse. Thus, we decided to espkue the set of values: 0. L 0.2 . .5, 
1.0 and 2.0 fmo. A size of idr-O.l um is representative of an unperturbed atmosphere, whereas rd-~2 
|im would mean that the entire taminator is rich in fi:esh ash. In going &aui r^^ =0.1 to 2 pm, the 
pihase EmctiDU in the forward diiection (calculated from klie tlieory for a unimodal log-Donnal size 
distribution, a geometric standard deviation of 2 and a lefiactive index of 1.41 M Ivlishchenko, 
w'ww.gi&sjiaEa.goV'ilaff''mrm&hclieiiko''t_n]atrix.litinl) augoients by -20. The phase fijnctioiis are 
calculated at a few wavelengths and interpolated in between. 
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Figure mentioned in the main text. 
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